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Abstract 

We study the general phenomenology of neutrinoless double beta decay in seesaw models. In 
particular, we focus on the dependence of the neutrinoless double beta decay rate on the mass of 

m 

^ \ the extra states introduced to account for the Majorana masses of light neutrinos. For this purpose, 

'x" ■ we compute the nuclear matrix elements as functions of the mass of the mediating fermions and 

; . estimate the associated uncertainties. We then discuss what can be inferred on the seesaw model 

^D \ parameters in the different mass regimes and clarify how the contribution of the light neutrinos 



should always be taken into account when deriving bounds on the extra parameters. Conversely, the 
extra states can also have a significant impact, canceling the Standard Model neutrino contribution 
for masses lighter than the nuclear scale and leading to unobservable neutrinoless double beta decay 
amplitudes even if neutrinos are Majorana particles. In particular, the decay rate is reduced by 
at least six orders of magnitude for masses of the extra states below 1 MeV in absence of extra 
contributions. We also discuss how seesaw models could reconcile large rates of neutrinoless double 
beta decay with more stringent cosmological bounds on neutrino masses. 
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I. INTRODUCTION 



At present, neutrino oscillations, together with cosmological evidence for dark matter and 
dark energy, present the only evidence for physics beyond the Standard Model (SM). Thus, 
models accommodating neutrino masses become an important component in the search for 
new physics and it is therefore fundamental to seek experimental answers to questions such as 
whether neutrinos are Dirac or Majorana fermions. In this context, one of the most promising 



OuB(3 decay), in 
-|9| for reviews. 



types of experiments is that of searching for neutrinoless double beta decay 

which a peak at the endpoint of the beta radiation energy is searched for, see 

Since this process is lepton number violating, its observation would imply that neutrinos are 

Majorana fermions [lO|. Present bounds on Oz//3/3 deca y fo r different nuclei have been se t by 



the Heidel 
Solotvina 
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berg-Moscow (llj JGEX (l2J, CUORICINO (l3|, NEMO 



and DAMA 



1^ 



y, ELEGANT |l6|. 



collaborations, as well as geochemical measurements [19|]. In 
the next few years, experiments like GERDA [20], EXO [21], SNO+ ^ and CUORE [231 



will search for a Oi//3/3 decay signal with unprecedented sensitivity. Furthermore, there are 



several other experiments proposed for the upcoming future, which include CANDLES 24 1, 



Super-N 
MOON 



MO 



25(, MAJORANA |26|, NEXT [27;], CARVEL [28;], COBRA |29|, DCBA 



3lj andXMASS |32|. 



m, 



Apart from the ongoing searches, there is a claim for a Oz//3/3 decay signal from (part of) 



the Heidelberg- Moscow collaboration 33|, |3J]. However, this measurement is in apparent 



conflict with other bounds on neutrino masses, in particular those coming from the influence 



of neutrino masses in cosmology 



35ll36|. 



Among the more popular models for neutrino masses, we find the different types of seesaw 
mechanisms. In the type-I seesaw 37H40|. the SM is extended by the addition of extra 
fermion singlets. The Majorana masses of such singlets do not violate gauge invariance and 
are presently not constrained. However, a very heavy mass is usually assumed and a Yukawa 
coupling between the singlets and the SM neutrino fields of the same order of magnitude 
as the SM Yukawas is introduced. Such a scenario provides a natural realization of small 
neutrino masses even if other possibilities, e.g., approximate lepton number conservation 



41 



42[, exist. 



In the current hterature on Oz//3/3 decay, it is common to study the effects of the SM 



neutrinos^ or the extra states introduced to account for their masses independently (see, e.g., 



Refs. 43N47l|). However, the interplay between both contributions, when combined, displays 
interesting phenomenology such as cancellations in certain regimes that is otherwise lost. 
By considering these contributions independent, the presence of some extra component to 
the Majorana neutrino mass different from the extra states and some degree of fine tuning 
is implicitly assumed. 

In the present work, we will compute the nuclear matrix element (NME) involved in 
the Oz//3/3 decay rate without any assumption on the neutrino mass mediating the process, 
describing in detail the approximations involved. The uncertainties associated to each of 
these approximations are also discussed, so as to estimate the total uncertainty on the final 



NME. The results of this computation are publicly available at Ref. 48|]. We will also discuss 
in detail the interplay between the contributions of the SM neutrinos and the extra states in 
seesaw models and discuss under which conditions they can be considered independent. In 
doing so we will cover the full parameter space and deduce what implications can actually 
be inferred on the models from observations. Finally, we will also comment on how the 
contribution of extra states with different mass scales could allow for a large Oi//3/3 decay 
rate even in presence of more stringent bounds on neutrino masses, such as those derived 
from cosmology. In particular, it will be shown that the Heidelberg-Moscow claim requires 
a tuning of only about 50 % amongst the extra contributions in order to be compatible with 
the present cosmology bounds. 

The remainder of this work is organized as follows: In Sec. [TTl we review the physics and 
assumptions used to compute the NMEs and perform these calculations in order to get their 
values as functions of the mass of the exchanged fermions. Next, in Sec. IIIIl we discuss the 
general phenomenology of the Oz//3/3 process and the approximations that are usually made 
in Oz//3/3 decay analyses. The realizations of the Ou/SP decay signal in the different types 
of seesaw models are treated in Sec. [IV] before we summarize and give our conclusions in 
SeclVl 



^ By SM neutrinos, we here mean the mass states which are predominantly composed of the SM flavor 
fields. Although these states can contain some admixture of the extra states introduced to account for 
neutrino masses, we will use this nomenclature for simplicity also in the remainder of the paper. 



II. COMPUTATION OF NUCLEAR MATRIX ELEMENTS 
A. Obtaining the rate of 0vf3f3 decay 

The Oz//3/3 decay is mediated by the weak Hamiltonian: 

Hw = ^{jL,J^')+^.c.. (1) 

where j^^ is the leptonic current, which consists of the electron coupled to left handed 
electron neutrino. The neutrino can be written as a linear combination of the light and 
heavy mass eigenstates Vi given by the mixing matrix U: 

3L^l = 67^ (1 - 75) ^eL, ^eL = ^ ^ei^^iL- (2) 

i 

On the other hand, the hadronic (nuclear) current J'^ can be obtained phenomenologically 
by imposing symmetry requirements to the more general combination that can be built with 
the available Lorentz vectors p^, p!^ and 7^, the neutron and proton four-momenta and the 
spin matrices, respectively. We also need to assume the impulse approximation, i.e., that 
nucleons in nuclei can be treated as free when dealing with the weak interaction. Then 
imposing Lorentz, parity and time-reversal invariance the nuclear current is given by^ 



jf = ^r' 



a 



flU 



9v{p h^ + Wm{p )- P,y - 9a{p )Yl5 - 9p{p )P^l5 

ZlTiN 



Vi/, (3) 



where p^ = Pn^Pp i^ ^^e transferred momentum from hadrons to leptons, rriN is the nucleon 
mass, vp represents a nucleon field and r~ is the isospin lowering operator, i.e., it turns a 
neutron into a proton. 

The form factors, gv, Qm, Qa and gp, are real functions of the Lorentz scalar p^. Their 
values at zero-momentum transfer are known as the vector, magnetic, axial and pseudoscalar 
coupling constants, respectively. Note that in single /3 and two-neutrino /3/3 decays only the 
vector and axial terms are usually considered, due to the small transferred momenta (< 1 
MeV). The magnetic and pseudoscalar couplings can be written in terms of the vector and 
axial ones by assuming the conserved vector current (CVC) and the partially conserved axial 



^ Second-class currents gs{p'^)p^ and gT{p'^)-§;^Pul5, for which there is no experimental evidence, will be 
ignored. 



current (PCAC) hypotheses [49|]. The CVC hypothesis also imphes that gv{0) = 1 in the 
nuclear medium. 

We now take the non relativistic approximation to the hadronic current. If terms are kept 
up to |p|/mAr (IpI ~ 100 MeV as will be discussed below) and the energy transfer between 
nucleons is neglected {E ~ p'^/2m]^), we are left with^ 

^f (X) = f2 ^n [a^'Jnip') + g'^Jnip')] ^(x - r.), (4) 



Ji=l 



where 

J'nip') = gvip'), 

Jnip') = ^9M{p')^f^ + 9A{p')cTr, - ^^(p2) P(P-^n) _ ^^^ 

Hence, we have a sum over all A nucleons of the nucleus, whose coordinates are denoted by 
r„. Note that nucleon operators present in the nucleon fields from now on will be included 
in the nuclear wavefunctions. 

The momentum dependence of the couplings is usually parametrized by the standard 



dipolar form [50(|, and takes into account that nucleons are not point particles but finite size 
bodies, i.e., the nucleon structure. Since in OufSP decay the neutrino is being exchanged in 
t-channel and the outgoing electrons have essentially the same energy, the energy exchange 
can be neglected and thus p^ ~ — P^- Then, the form factors look like 

9v{p^) = —^ -2, 9m{p^) = (^p - /in) 9v{p'^), 

(i-^ft) 



(-0 



where 771^^ is the pion mass and fip and /i„ denote the proton and neutron anomalous magnetic 
moments, respectively. The values of the cutoffs of the vector and axial nucleon form factors. 



Ay = 0.85 GeV and A^ = 1.09 GeV, are taken from experimental observations 5l|,l5^- Their 



effect is to weaken the couplings for large momentum transfers, i.e., they avoid contributions 



■^ Nuclear recoil terms also come at first order in l/rriN, being proportional to Pp + p„ instead of p. 
However, their leading contribution is suppressed an extra order of magnitude because of their odd-parity 
character, which requires electron p-waves. Hence these terms will be neglected. Their contribution will 
be discussed when referring to the electron s-wave approximation. 



arising from nucleons being too close to one another. These form factors are commonly 
denoted as the finite nuclear size (FNS) terms. 

With this Hamiltonian the rate of the Oz//3/3 decay can be calculated by using the second 
order Fermi's Golden Rule 2|: 

dTo^fip = 27r ^ \Roui3pf S{ei + £2 + Ej - Ei)dVte^dVte^, (7) 

spin 

where the transition amplitude is given by 

2 



Roupp = f -^ j j d^ j dy-^ (1 - P12) X 



>p [Nf] ei, 62! J'fi{^)]L^,i?^) \Na] ei, z/j) {Ng, d, Uj\ Jl\y)3Lpiy) \Ni) 

^ Uj + Ea- {E, + £1) ■ ^ ^ 

The operator P12 is included to fulfill antisymmetry for the electrons, whose energies are 
denoted by £1, €2- The energy of the virtual neutrino Vj is denoted by Wj, and Ea is the 
energy of the virtual intermediate nuclear state liV^). 

The leptonic part of the numerator in Eq. (jHj) can be written more explicitly as 



-I Y}JlA^)l, (1 - 75) / '.^ o [t + va,) C^ (1 - 75) l,e^{y) 



3 

74„ pip-{x-y) 



= -'1^ ^ej^j J -^^-^2eix)^^ (1 - 75) 7pe^(y) 

= - ^ E ^ejm, I -^^-^-^e{x)^, (1 - 75) 7pe^(l/). (9) 

^ J (27r) Wj 

This term turns out to be proportional to the neutrino masses ttij because of the left handed 
character of both leptonic currents. The transition amplitude is then 

i?o./3;3 = ^^^^ I rfx I rfy (1 - P12) e(5i, x)7, (1 - 75) 7pe^(^2, y) X 

Y "^e 3\ (0) J 27r2 Wj ^ OJj+^la-l {si -62) 

The axial coupling (7^(0), the electron mass irtf. and the nuclear radius R have been intro- 
duced for convenience and to make the second line in Eq. flTOj) dimensionless. For the nuclear 
radius we have taken R = 1.2 A^^^ fm. Using energy conservation, we have also rewritten 
the denominator using the new parameter 

^^^^ET-^iE^ + Ej), (11) 



which gives the relative energy of the (virtual) state of the intermediate nucleus with respect 
to the mean energy of the initial and final states. Typical values of this parameter for the 
different decays are ~ 10 MeV [ij. 

In the following, two approximations will be made: 

• Closure approximation. 

• Approximation of 0"*" final states and electrons emitted in s-wave. 

The first of these takes advantage of the high momentum of the virtual neutrino |p| ~ 
100 MeV. This nuclear scale comes from the integral over the transferred momentum of 
Eq. ( ITOj) . For light neutrinos, the integrand of this equation is approximately proportional 
to IpI /(IpI + fia) in radial coordinates. Hence, momenta < 10 MeV will be disfavoured in 
the transition. In addition, if we recall the form of the FNS terms in Eq. (|6]), we see that 
transferred momenta above the cutoffs ~ 1 GeV will be suppressed as well. On the other 
hand, the leading contribution of the exponential term will arise when p- (x — y) c^ 1. Since 
nucleons are typically few fermis (femtometers) apart in nuclei, the nuclear wavefunctions 
will select the preferred momentum for the virtual neutrino to be |p| ~ 100 MeV. Thus, we 
would expect it to be the typical virtual neutrino momentum of the decay. This result has 



been confirmed by explicit calculation 53l. |54|. 



In the case of heavy neutrinos, the transition operator will now have stronger preference 
for larger momenta |p| ^ nij, since Uj = \/fn'j + p^ appears in the denominator in Eq. ( ITO|l . 
Therefore, the tendency of the nuclear interaction for |p| ~ 100 MeV can be overcome result- 
ing in large transferred momenta that would imply internucleonic distances much shorter 
than ~ 0.1 fm. Again, such a distance is very suppressed by taking the FNS terms into 
account, i.e., the nucleon structure information. Hence, in these cases, the reduction due to 
the FNS effects is very large so that we end up with a |p| value of a few hundreds of MeV at 
most, which is also the expected value for a process taking place between nucleons in nuclei. 

In any case, the term {ei — 62), which can amount up to a couple of MeV and vanishes 
on average, can be safely neglected. Moreover, the intermediate state energies E^, which 
can differ from one another by a few MeV, can also be replaced by an average value (£'"*). 
Thus, only a common parameter 

/x, ^ /i ^ (i?-) -liE, + Ef) (12) 



is required [l| . With the removal of the dependence on the actual energy of the intermediate 
states, it follows that they are no longer needed in the calculation, since the closure relation 
can be applied. Thus, the hadronic part in Eq. (TTUj) is now 

E V W , ^ 77—7 E ^/ -^Li^) Na) {Na Jf y iV. 

= ^— (iV/l^{:^(x)J?(y)|iV.). (13) 

This closure approximation has been shown to be correct to more than 90 %, using the 



quasiparticle random phase approximation method 55|, to be presented in Sec. IIIBI 

As for the limitation of our study to transitions to 0+ final states, and to cases where 
electrons are emitted in s-wave, corrections are expected to be of the order of 1 % at most. 
In the case of p- waves, they are suppressed to s- waves by an order of magnitude at least |2|. 
Moreover, since they have odd-parity, they need odd-parity terms in the current to couple 
to a J"*" final state and these only appear at (|p| /rnN), see note on Eq. (jlj). Thus, in the 
end we have a contribution of the order of 1 % of the leading ones. On the other hand, since 
all the final nuclei of Oz//3/3 decay have even number of both protons and neutrons, their 
ground states are always 0^. Any transition to an excited final state will be suppressed by 
a phase space factor, which in this case is approximately proportional to Q^^, where Q/3^ 
is the energy available for the decay. Considering this factor, the only other low-lying final 
states of interest are 2+ excited states. Moreover, apart from the phase space suppression, 
these final states also need electron p-waves due to angular momentum coupling, since two 
electron s-waves can only couple to angular momentum or 1. Hence, these transitions can 
also be safely neglected. 

Within these approximations, the transition amplitude can be written as 

V t/2.^ (0,1 Vr-r-—^ /iPe^p(rn-r™)_^Wp!L|o.\ (U) 

J n,m ^ ^ ' J \ J I / 

where the index contractions and the x and y integrations have been performed. The 
operator Qnm (p^) is the result of the product of the nuclear currents given in Eqs. (jlj) and 
(j5j), and has Fermi, Gamow-Teller and Tensor contributions: 

nnmip') = -h^ip') + /^^'^(P')^n " <Tm ' k"" {p')S^^^, (15) 



with the tensor operator SJJ^ = 3 (p ■ cr„)(p ■ cr^) — cr„ ■ cr^. The functions /i(p^) can 
be labeled according to the terms of the hadronic current [see Eq. ([3])] from which they 
originate: 

h^ip') = hUp') + hUp') + hlUp')- (16) 



Their explicit form can be found in Ref . 46| . It was in this work that the importance of the 
non leading terms [i.e., all but h^^ and hyy in Eq. (fT6|) ] was first shown. They are referred 
to as higher order components of the nuclear current (HOC). Notice that, since these terms 
are of orders |p|/m7v and {\p\/mN) in the current,^ their contribution will be enhanced for 
larger transferred momentum, i.e., for heavy neutrinos. 

It is easy to insert Ro^pp into Eq. ([7]) and derive an expression for the Oz//3/3 decay rate: 

(17) 



2 



tJOl 



In 2 



[rrij^ 



Here, Gqi is a well known kinematic factor, and comes essentially from the leptonic degrees 
of freedom. It can be written explicitly as 

Goi = '^4^5^°^i^7' / ^°(^' ^i)^o(^' ^2)gig2ei52<5(£i + £2 + ^/ - E,)de^de2d (qi ■ qs) , 

(18) 

where F^iZ^e) are the so-called Fermi functions, with Z the proton number and qi the 
electron momenta. The Fermi functions also depend on the nuclear radius R and their 
explicit form is 

2(71-1) ip/ ,-N|2^^y. ^,. - .[7~7Z^2 "^^ 



Fo(Z,5) = p^^^^^^) {2qKY^'^^-'^\V{n,^%yte-y- ^, = ^j\-{aZ)\ y ^ . 

\l9) 
The quantity M^'^^^{inj) is the nuclear matrix element, which takes into account the initial 
and final nuclear wavefunctions and the transition operator. This operator originates from 



The terms labeled hpp and Hmm come at order (|p|/7nAr) . The reason to keep these second order terms 
is the enhancement of the coupling constants gp and gn due to the factors 2mN |p| / (p^ + w^) — 7 
and (/Zp — /x„) = 3.70, respectively, see Eq. ([6]). In each of these terms the corresponding factor appears 
squared. 

9 



both the nuclear currents and the virtual neutrino. The NME is given by 

M»-(,„,) ^ (0;| E r,r, JL. I * ,.p.<..-...^iW£^ |0.) . (20) 

When the integral over p is performed, we obtain 

M'^^^im,) = (0+1 5^ r-r„ (-V^(r) + V^^(r)a„^^ - V^{r)S:^) |0+> , (21) 

n,rn 

where r = |r„ — r^l is the distance between the decaying neutrons and the V{r) are the 
so-called neutrino potentials. Before the radial integration over |p| they are given by 



^.W = ---f77^/ J2(|p|r) ;7 p^c^IpI, (22) 



where j„(a;) are the spherical Bessel functions. 



B. Calculation of the nuclear wavefunctions within the Interacting Shell Model 



As can be seen from Eq. (12T]) . a key ingredient in the calculation of the NMEs are the 
wavefunctions of the initial and final nuclei. This is a complicated nuclear structure problem 
which cannot be solved in the complete space, i.e., taking into account all neutrons and 
protons of the corresponding nucleus in all their possible configurations. Thus, truncated 
valence spaces and effective interactions are used to solve the nuclear many body problem. 
As a consequence of this, a fully consistent treatment would demand regularizing the Ou0B 



decay operator in Eq. ( 12T]) using the same prescription as for the bare nuclear interaction 56 1. 



This has only been performed very recently 57, 



Instead, we will simplify the problem keeping the bare Oz//3/3 decay operator and including 
new correlations [called short range correlations (SRC)] in the calculation via a general 



prescription. The findings of Refs. [57| and [58| show that for light neutrinos, the effect 



of these correlations is rather moderate once FNS terms have been taken into account, of 
the order of 5% correction to the NME. In order to implement these SRC one needs to 
assume some prescription and the most commonly used are either a Jastrow-type function 



59| or a unitary correlation operator method (UCOM) transformation [60]. The NME is 



thus transformed as 

(0+1 O'"^" \Ot)sRc = (0/1 U^O'-'^^U |0+> = (0+1 0°'^'^^ |0+> , (23) 

10 



where U is either a Jastrow-type function or a UCOM transformation. The actual 
parametrizations can be found in Refs. 58|] and 6l|, respectively. In our calculations have 
used the UCOM prescription, even though similar results are expected within the Jastrow 
approach. Note that these SRC terms, which do not have much importance for light neu- 
trinos, will be more relevant for heavy ones, which require shorter distances between the 
decaying nucleons. 

Mainly two different methods are used to obtain the NMEs for the Oz//3/3 decay, the 
quasiparticle random phase approximation (QRPA) [J, |5| and the interacting shell model 
(ISM) 62|, |63| . The QRPA includes relatively large valence spaces but is not able to comprise 
all the possible configurations. On the other hand, the ISM is limited to smaller configuration 
spaces, but all possible correlations within the space can be included. 

When comparing the NMEs obtained by both methods a clear disagreement is found, 
the ISM values being about 1.5-2 times smaller than QRPA ones 54|. This disagreement is 
not overcome when taking into account the estimated errors of both calculations and it is 
under discussion whether the difference is mainly due to the lack of correlations of QRPA 
calculations (which reduce the value of the NME), the small valence space used in the ISM 



(to be discussed in Section III CI) or both [53 



631. However, it must be stressed that due to 



the theoretical effort made over the last years this disagreement is now much less severe than 
it was five years ago, and studies using recently available experimental information for the 
decay of ^^Ge suggest that the present situation can be improved 



M, 



651]. This applies to 



NMEs obtained with light neutrino exchange. In the heavy neutrino case, since the available 
QRPA results are rather outdated |i46|, it is probably not meaningful to compare them with 
those of the present work. However, the same relation between ISM and QRPA results of 
the light neutrino case is to be expected, since the difference between these methods lies on 
the calculation of the wavef unctions, and therefore it should not be very much dependent on 
changes on the transition operator, for whom both methods give an equivalent description 



54| 



In this work, we have used ISM nuclear wavef unctions. Following the considerations of 
Sec. Ill At we have performed calculations for the tjij dependent Oz//3/3 decay NMEs of the 
emitters ^^Ca, ^^Ge, ^^Se, ^^^Sn, ^^^e and ^^exe, using the ISM coupled code NATHAN [ej, 
ideally adapted for the calculation of 0"*" states. Full diagonalizations are accomplished within 
different valence spaces and effective interactions. For instance, the decay of '^^Ca is studied 



11 
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FIG. 1: The nuclear matrix element dependence on the neutrino mass. The data depicted in this 
figure is available at Ref. [48 1. 



in the pf major shell, where the KB3 interaction [66| is employed. For the case of ^^Ge and 
^^Se, the valence space consisting on the lp3/2, O/5/2, lpi/2 and Ogg/2 orbits is diagonahzed 
using the GCN28.50 interaction [S^l- Finally the OgT/2, l<i3/2, l'^5/2! 2si/2 and O/111/2 valence 
space and the GCN50.82 interaction JSJ] are used in the decays of ^^^Sn, ^^°Te and ^^^Xe. 
The results of the computation of the NME as a function of the neutrino mass are shown 
in Fig. [H which is in agreement with the findings of Refs. [67] and 43|. As can be seen 
from this figure, the dependence of the matrix element on the nuclei is mild and all the 
curves show a similar behaviour. We will discuss qualitatively the observed dependence of 
the NME on the neutrino mass as well as the phenomenology associated to the different 
mass regimes in the next sections. 



C. Estimate of the uncertainties of the NME within the ISM 

As we have pointed out several times throughout the present section, since the com- 
putation of the NMEs can be rather challenging, a number of approximations need to be 
made. In consequence, the results obtained will have some uncertainties, which we will now 
estimate. 



12 



First of all, when dealing with the ISM, we have to worry about the valence space and 
the effective interaction used to obtain the nuclear wavefunctions. The effect of having 



larger valence spaces was analyzed in Ref. 68|, with the result that the NMEs increased. 



within a conservative analysis, by ~ 15 — 20 %. The same number was obtained in a QRPA 
calculation when it was quantified the effect of the orbits absent in a ISM valence space 



6J]. Of the three valence spaces employed in the present work (see Sec. IIIBp . only that 



corresponding to the ^^Ca decay was not explored in Ref. 68|, so that we can consider the 
former ~ 15 — 20 % increase as a general estimate of the uncertainty due to the configuration 
space. As for the nuclear interaction, a moderate dependence of ~ 5 — 10 % was found both 



in Refs. 68| and [65|. Since they study nuclei in different regions and we can consider all the 
effective interactions employed of similar accuracy, we will also take this figure as general. 

Until recently, another considerable source of uncertainty were the SRC. However, it 
now seems that their contribution is rather small and that proper UCOM or Jastrow-t ype 



parametrizations can take these terms very well into account, with a precision of ~ 5 % 57, 



For its part, the variation in the NME due to different but reasonable values of the 



cutoffs appearing in the FNS terms [49|, |69| is also very small, less than 5 vo. 

In addition, due to the fact that the NMEs have been calculated using the closure ap- 
proximation, we have to include an additional error of 5 — 10 % to our results, as suggested 
by QRPA calculations [55!] • This is in agreement with the very soft variation that is seen 
in the NME as the parameter ^ is modified. Furthermore, we quantify the possible effect 
of missing terms in the operator (next order terms in the current, nucleon recoil, p-wave 
emitted electrons) by an additional uncertainty of less than 5 % in our results 



7d, 



. t is currently under discussion whether the axial coupling should be quenched or not 



65 



71| . In the ISM calculations presented here we take (7^(0) = 1.25, i.e., we do not quench 
it, contrary to what is required by the single /3 and two- neutrino (3(3 decays, where a pure 
Gamow- Teller operator appears and its value has to be quenched to 5'a(0) = 1.00. However, 
in the Oz//3/3 decay case the operator is more involved due to the extra radial dependence 
introduced by the virtual neutrino. Moreover, in the case of the pure Gamow- Teller J^ = 1"*" 
channel, it is not dominant in the Ou^P process, and depending on its relative sign, quenching 
it may result even in an enhancement of the NME. Until this issue is explored in more detail, 
we will take the most conservative option, allowing for a full quenching of the axial coupling 
and also for the quenching only of the J^ = 1+ channel. Under these assumptions, taking 

13 



into account that the Fermi part of the NME accounts for 10-15% of the full NME in our 
ISM calculations and is never quenched, we estimate a ~^3q % error due to this effect. 

Notice that the uncertainty in the valence space only moves the estimate up, the effect of 
axial quenching essentially moves it down, while the remaining contributions are expected 
to be Gaussian-distributed. Even though some of these errors may be correlated in a rather 
complicated way, as a first approximation we will take them as independent. Altogether, 
adding every contribution in quadrature we expect an overall uncertainty in the final NME 

^r +25 % 
Oi ~-35 %• 

The above analysis applies to the case of light neutrino exchange. For heavy neutrinos, 

the NMEs get very dependent on the SRC and FNS treatments. In this case, a similar study 

to the one above gives a 15-20% for SRC uncertainties and 10% for FNS ones, which would 

lead to a final ~!!^4q % uncertainty. However, one should take this number cautiously until 

ihe accuracy of the FNS approach for such heavy exchanged particles is firmly established 



72|. 



As an example, we will consider the case of the ^^Ge decay. In Ref. 65|] the NME 
was obtained with different effective interactions and SRC, obtaining the interval 2.81 < 
]\^ov/3/3^g^ < 3.52. If we take into account the further uncertainties of the valence space, 
the FNS, the closure approximation, the next order hadronic current terms and the gA{0) 
quenching, we end up with 2.11 < M^'^^^{0) < 3.98. This result will be used in Sec. IIVI to 
derive bounds on the neutrino masses from the Oi//3/3 decay process. 

III. GENERAL PHENOMENOLOGY 



According to Eq. fITTj) . the contribution of a single neutrino to the amplitude of 0^(3/3 
decay is given by 

Aocm,f/f,M°'^^'^(m,), (24) 

where rrii is the mass of the propagating neutrino and M^'^^^inii) is the nuclear matrix 
element that characterizes the process and depends on the nucleus that undergoes the Oi//3/3 
transition. Figured] shows two distinct regions where the behaviour of the NME as a function 
of the neutrino mass changes from almost constant up to rrii ~ 100 MeV to decreasing 
quadratically as the neutrino mass increases beyond 100 MeV. This behaviour is easily 
understood: the neutrino can be characterized as light if mf <^ |p^| or heavy if mf ^ \p^\, 
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which would mean that the neutrino propagator in the NME would be dominated by p^ or 
mf, respectively, where p is the momentum exchanged in the process. As already mentioned 
in Sec. IIIAl in the Qv(3(3 decay p^ ~ — p^ ~ —(100 MeV)^. We will therefore define two 
regimes: 

• The light neutrino regime: For irii < 100 MeV, where the neutrino propagator is 

1 



p^ — mf 



1 m2 


(m^ 


p2 pA 


\f 



(25) 



and hence, the NME is maximum in this regime and is almost independent of the 
neutrino mass: 



M°''^^(mi) = M°'^'^''(0) 



i + ^ + oi^ 

pZ y p4 



(26) 



• The heavy neutrino regime: For rrii > 100 MeV where the NME decreases as 
m^'^ providing an extra suppression to its contribution to the 0^(3(3 decay amplitude 
because of the neutrino propagator: 

^^ = -A + ^f4V (27) 

p'^ — m^ mf \'^^i J 

In principle, one could expect a more involved behaviour or even a resonance if p^ ~ mf. 
However, since the Oi//3/3 transition does not occur through an s-channel type diagram, the 
characteristic momentum transfer has p^ < with |p^| ^ (100 MeV)^. Thus, the nuclear 
matrix element does not exhibit a resonant mass and the transition between the light and 
heavy neutrino mass regimes is relatively smooth (c./.. Fig. [1]). We will therefore not define 
a third transition region between the heavy and light regimes, as is sometimes done in the 
literature, since there is no new phenomenology associated to it. 

The usual bound derived from the Oz//3/3 process in the literature is obtained summing 
over the active neutrinos, implicitly neglecting the contribution of extra degrees of freedom. 
With this assumption the only contribution comes from the neutrinos in the light regime, 
with iTLi <^ p^, which results in 

3 3 3 ._2^ 



Ao.,, = E ^^ °^ E ^^UlM'^^^im.:) = M°^^^(0) ^ m^U^ + O (^) , (28) 

where Yli''^i^ei is the well-known expression used for the "effective Oi//3/3 decay neutrino 



mass" : 



mfjl3 = micl^cl^ + m2sl2cl^e^'°'^ + m:>,sl^e^^°"' , (29) 
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FIG. 2: Allowed regions by neutrino oscillation data at 2a for m^/j as a function of the mass 
of the lightest neutrino mi, for normal (blue) and inverted (red) neutrino mass hierarchies. The 
cosmologically disfavoured neutrino masses from cosmology are also shown (green band) as well as 
the Heidelberg-Moscow claim for Oz^/3/3 decay (yellow band). 



cos Uij , Sjjf 



sin 6ij, 6ij are 



where rrii are the masses of the neutrino mass eigenstates, Cij 
the neutrino mixing angles and ctj are combinations of the Majorana and Dirac phases. It 
is important to note that this expression holds only when the SM neutrinos dominate the 
OufSP process. However, if the SM is not extended, the Majorana mass required for the 
Oz//3/3 transition is forbidden. As we will discuss, all extensions of the SM that induce a 
Majorana mass for the SM neutrinos imply the inclusion of extra degrees of freedom that 
can contribute to the OuPP process and should be added to the SM decay amplitude. 

Under the assumption that only the light neutrinos give a significant contribution to this 
process. Eg. ( l29l) can be combined with the present constraints on neutrino masses and 
mixings [73| to derive Fig. |2l where the allowed value of m^/3 as a function of the mass of 
the lightest neutrino mi is shown. However, as will be discussed in Sec. llVt this assumption 
is not always valid and Fig. |2] does not always provide an accurate description of the OuPP 
decay. In the same figure we also show the Heidelberg-Moscow claim for Oz//3/3 decay and 
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the bound on the mass of the hghtest neutrino from cosmology arguments |35l . |36| . Notice 
that the bounds from cosmology apply to the SM active neutrinos only. For extra sterile 
neutrinos, the cosmology bounds would depend on their mixing with the active ones, which 
would determine their abundance. Indeed, sterile neutrinos with masses ~ 1 keV larger 
than the cosmology bound are actually considered good candidates for warm dark matter 
74J-|76| . As can be seen from the figure, there is a tension between the Heidelberg-Moscow 
claim, the contribution to Oz//3/3 from SM neutrinos and the present bounds on their mass 
from cosmology. In Sees. IIVA3I and IIVDI we will discuss possible solutions to this tension. 

IV. APPLICATIONS TO SPECIFIC MODELS OF NEUTRINO MASSES 

In this section we analyze the contributions to Oz//3/3 decay of the different mechanisms 
that lead to Majorana neutrino masses and can therefore induce the required lepton number 
violation for the process. We will discuss here the tree-level realizations of the Weinberg 
d = 5 effective operator: 

^ f-n A.A (l\ 



^ y^u.^ J ^, L^j+H.c. . (30) 

Here, (f) denotes the SM Higgs field, which breaks the electroweak (EW) symmetry after 
acquiring its vacuum expectation value (vev) f , A is the scale of new physics that gives 
rise to the operator and we have used the definition cj) = zr20*. This is the only d = 5 
operator that can be built from the SM particle content respecting both gauge and Lorentz 



invariance [77|]. Since the low-energy effects of physics beyond the SM can be encoded in an 
expansion of effective operators of d > 4 suppressed by inverse powers of A*^"^, it is a natural 
expectation that this sole d = 5 operator will be the least suppressed one. It is then very 
suggestive that one of the few evidences we have for physics beyond the SM is the existence of 
small, but non vanishing, neutrino masses. Indeed, after the Higgs field develops its vev, the 
operator of Eq. ( !30|) induces a Majorana mass term for the SM neutrinos Ca/3(f ^/A)z/^^z//3i, 
suppressed by the scale A. 

There are three different extensions of the SM particle content that lead to the operator 
of Eq. fl30l) after the extra mediators have been integrated out. They are known as "seesaw" 



mechanisms of type-I 



37 



40|], where the heavy particles are fermion singlets, type-II [781 - 

, where the SM is extended 



82|, where scalars triplets are included, and type-Ill [83 
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FIG. 3: Feynman diagrams contributing to the Oi^/3/3 transition rate from the exchange of fermions 
(left) and scalar triplets (right). 

by fermion triplets. All these extra degrees of freedom, required to induce the Majorana 
nature of the SM neutrinos, can also contribute to the OufSP process. The contribution of 
neutral fermions, such as the singlets and triplets added in the type-I and III seesaws, to the 
OuPP decay rate is depicted on the left side of Fig. [Sj in particular this includes that of the 
light active neutrinos. The contribution of the scalar triplet of the type-II seesaw (see, e.g., 
Ref. 87|) is depicted on the right side of Fig. |3l In both diagrams the W lines can also be 
exchanged for the physical singly-charged scalar present in the type-II seesaw. In principle, 
the contributions of the light active neutrinos and those of the extra degrees of freedom that 
are introduced should be combined. This is especially so in the case of the type-I seesaw, 
since important cancellations 



90| are present in certain regimes [9l|| that are missed 



if the constraints are placed separately. If they are assumed to be independent, an extra 
contribution to the neutrino mass beyond those extra states is implicitly assumed. Moreover, 
taking into account the relations between the high- and low-energy parameters allows the 
derivation of stronger bounds on the former through the active neutrino contribution. It 
is then important not to neglect it, since the naive constraints stemming directly from the 
contribution of the extra degrees of freedom are generally much weaker. 

A. Type-I seesaw models 

In this section we will discuss the phenomenology of Oz//3/3 decay when extending the 
Standard Model with fermion gauge singlets, i.e., right handed neutrinos z/^j. The Standard 
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Model Lagrangian is then extended as^ 

=^ = ^SM - -17-{MN)ijy% - {Yr,)iaT^^^L^ + H.c. . (31) 

After the Higgs develops its vev the neutrino mass matrix is 

f Ynv/V2\ , , 

M.= \ ^ ^ . (32) 

This mass matrix can be diagonalized by a unitary mixing matrix IJ: 

f/* diag {mi,m2, ...,m„} ?7^ = M,^ . (33) 

We therefore have n neutrino mass eigenstates with masses mj and mixings f/gj with the 
electron. Out of these, at least three mass eigenstates must be very light and form the main 
components of the active neutrinos, whose number is measured by the invisible decay width 



of the Z 92] • On the other hand, the masses of the extra states are not determined. We thus 
have two contributions to the amplitude of Oi//3/3 decay, one from the light active neutrinos 
and another from the extra degrees of freedom: 

A ex ^ mif/2.M°'^^^(mi) + ^ mj\Jl,M'''^\mj), (34) 

i I 

where we have used capital letters to denote the mass index of the mostly sterile states 
and lowercase letters for that of the mostly active states. Depending on whether the extra 
mass eigenstates fall in the light or heavy neutrino mass regimes we can further split their 
respective contributions to the amplitude: 

light light heavy 

A oc ^ mi?7f,M°"^^(mi) + ^ m,?7f,M°'^^^(m7) + ^ m,t/2,M°"^^(m/). (35) 

i I I 

We can now distinguish three cases exhibiting very different phenomenologies depending on 
the mass regime of the extra mass eigenstates: 

1. All extra mass states in the light regime 

In this scenario all the mass eigenstates are lighter than 100 MeV. In principle, this does 
not allow to explain the smallness of neutrino masses through the naive seesaw mechanism. 



^ Here, as in the following, we do not write out the kinetic terms of the new fields, which are assumed to 
be of the canonical form. 
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FIG. 4: Nuclear matrix element cancellation in the light mass regime. The example is taken for 
the ^^Ca decay, but the same quadratic dependence is seen for the remaining nuclei studied. The 
value of the parameter is a = 2.51 • 10^'^ MeV^^. 



i.e., with O (1) Yukawa couplings. However, since the value of the parameter Mtv in Eq. (l3Tll 
is not restricted and it is technically natural for it to be small (for vanishing Majorana mass 
term the B — L symmetry is recovered), we believe this is a possibility worth exploring even 
if less appealing than the canonical type-I seesaw scenario. 

Notice that, if all neutrinos belong the the light regime, Eq. f l33|) implies 



light light 

i I 



(36) 



since the left-left entry of the mass matrix in Eq. ( I32l) vanishes (it is forbidden by the SM 
gauge symmetry). Thus 



light 



light 



A oc ^miUlM'"'^^{mi) + ^miUl,M'"'^^{mi) 



light 



Y^^iUli (m0'^''''(0) - M°'^'^'3(m,)) . 



(37) 



Indeed, for neutrino masses in the light regime, the nuclear matrix elements are basically 
independent of the neutrino mass (see Fig. [ID M^^^i^inii) ~ M'^''^i^{mi) ~ M°''^^(0) and 
therefore the rate of 0^(3(3 decay in Eq. (|37|) is very suppressed when all mass eigenstates 

20 



10-3 

10-s 

10-12 
10-15 

10-*^ 




Complete y////\ 
Extra V only 



1 10-^ 

mi (MeV) 



FIG. 5: Bounds from CUORICINO on the extra neutrino mixing from Ou/3/3 decay in ^^Te, with 
a 90 % CL half-hfe [13(]. We have assumed the extra neutrinos are non hierarchical and show the 
bounds as a function of their common mass. We compare the case in which the contribution from 
the SM neutrinos is properly taken into account (striped area) to that in which only the extra 
contribution is considered (above the red line). 



are lighter than ~ 100 MeV. This works in a way similar to the GIM suppression in flavour 
violating processes [93]. Indeed, a process converting from flavour a to flavour /3 is GIM 
suppressed due to the unitarity relation J2i UaiU^i = 0, so that only the mass of the propa- 
gating particle does not make the cancellation exact and a suppression of Arn^/M^ is found, 
where Am^ is the mass squared difference between the different propagating particles. In 
our Oz//3/3 process we have Eq. (I36ll operating a similar cancellation and again it is the differ- 
ent neutrino masses in the NME of Eq. ( 1371) that would prevent a full cancellation leading 
to a suppression of Am^/p^ with |p^| ~ (100 MeV)^ [see the second term in Eq. fl26|) and, 
e.g., Ref. [QJ] ]. The Arn^/p^ dependence of the M°^''^(0) - M^''l^^{mi) term that drives the 
contribution to Oi//3/3 in Eq. fl37|) is depicted in Fig. HI As expected, the two contributions 
cancel up to a factor mj/p^ with \p'^\ ~ (100 MeV)^ and deviations from this behaviour start 
to be non negligible for mj > 1 MeV. 



We want to emphasize that phenomenological analyses (see, e.g., Refs. 43| and 44| ) 
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that use the non observation of OufSP decay to derive bounds on the mixing of an extra 
hght neutrino with mass around ~ 100 MeV or below neglecting the contribution of SM 
neutrinos, implicitly assume that the cancellation described in Eq. f l37p does not take place. 
This would be the case if some extra contribution to the neutrino masses is present and 
Eq. (136 p is consequently modified. Examples of this situation, which imply some degree 
of cancellation between the different contributions to neutrino masses, will be discussed in 
Sees. IIVA3I and IIVDI In the absence of these extra contributions, the remaining leading 
term, GIM suppressed as Am'^/p'^, can then be used to derive the corrected bound on the 
mixing of the extra state. In Fig. [S] we show this bound using the constraints on OuPP 
decay in ^^OTe from CUORICINO, with a 90 % CL half-life pj] and compare it to the one 
that would be obtained neglecting this cancellation and considering only the extra neutrino 
contribution. It is clearly seen that the bound indeed weakens when mi < 100 MeV. In 
particular, for mj < 1 MeV, the M^"^^ {0) — M^'^^'^ {mi) suppression in Fig. IHbecomes almost 
10^^, so that the 0^(3(3 decay rate would be suppressed by six orders of magnitude and for 
mi < 100 keV the bound on the mixing from Fig. |5] becomes meaningless since Yli ^e/ < 1- 
Clearly, if the extra states that accommodate neutrino masses are all in this regime the 
OuPP decay becomes experimentally inaccessible even if neutrinos are Majorana particles. 
An important sub-case is when M^ = with three extra neutrinos, which corresponds to 
the case of Dirac neutrinos. In this scenario, the cancellation is perfect since the left and 
right handed states are exactly degenerate, which corresponds to the fact that there is no 
lepton number violation for Dirac neutrinos. 

In order to derive the constraints of Fig. [5] we have assumed as an example that all the 
extra states are degenerate in mass (or that there is only one extra state). However, a similar 
GIM-like cancellation would also be present when assuming a different hierarchy as long as 
all the extra states are in the light regime. The extra contribution in the light regime can 
only dominate over the light active one and avoid the GIM-like cancellation if the model is 
extended with other heavier states. This situation also implies a certain degree of fine-tuning 
between the extra contributions and will be discussed in detail in Sees. IIVA3I and IIVDI 
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2. All extra mass states in the heavy regime 

This is the simplest situation and the one imphcitly assumed when using Eq. (129!) or Fig. [2] 
to forecast the rate of Oz//3/3 decay. It also corresponds to the canonical implementation of 
the type-I seesaw mechanism, where the Majorana mass of the right handed neutrinos in 
Eq. (pT]) is assumed to be very large Mjv ^ Y]^v /^/2 so that the smallness of the active 
neutrino masses is naturally explained and extra heavy eigenstates are obtained. Since the 
extra states would fall in the heavy neutrino mass regime their nuclear matrix elements are 
very suppressed: M°^^'^(m/) ^ M^'^'^^{mi) (see Fig. [1]). Furthermore, Eq. f p3|) implies that 

light heavy 

E^.f^e.+ E"^^f^e/ = 0, (38) 

i I 

SO that the Oi//3/3 decay amplitude is 

light heavy 

A oc ^mi[/2.M°"^^(mi)+ ^ mjU^jM^'^^imj) 

i I 

heavy 

I 

heavy light 

^ - ^ mjf/2,M0'^^^(0) = ^mi[/2.M0'^^^(0). (39) 

/ i 

The contribution from the light active neutrinos thus dominates and Fig. |2] provides an 
accurate prediction for the Oz//3/3 transition rate. 

Notice that, using the contribution of the extra states mjU^jM'^'^^^{mj) to the Oi//3/3 
process in order to derive a bound on their mixing Uei would lead to rather weak constraints 



since it is very subleading (see, e.g., Refs. |45l-l47j ) . On the other hand, Eq. (155]) can be used 



to express the dominant light neutrino contribution as a function of the heavy parameters: 
—miU^jM^'^^^{0). This expression then allows to derive a much stronger constraint on Uei 



95| . This constraint is also shown in Fig. [5] when mj > 100 MeV and compared to that 
when only the contribution of the heavy neutrinos is considered. In order to derive the 
constraints of Fig. |S]we have assumed as an example that all the extra states are degenerate 
in mass (or that there is only one extra state). However, a similar behaviour is found when 
a different hierarchy is considered as long as all the extra states are in the heavy regime. 
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3. Extra mass states in the light and heavy regimes 

In this scenario we would have the full contribution to Oi//3/3 decay of Eq. (135|) and the 
constraint from Eq. f l5^ : 

light light heavy 

Y, ^^Ul + £ miUl, + Y. niiUJ, = 0. (40) 

i I I 

As discussed above, for the neutrinos in the heavy regime, the NME receives an extra 
suppression to their contribution to the Oz//3/3 decay rate and the leading terms stem from 
the light states: 

light light heavy 

A oc ^ nXiUlM'^^^^imi) + ^ mj[/2,M°'^'^^(m,) - - J^ m,f/fjM°"^^(0). (41) 

i I I 

However, in this case the GIM-like cancellation is prevented since the heavy contribution is 
suppressed. This scenario thus offers the richest phenomenology. In particular, it is possible 
to satisfy Eq. (HOj) even in a situation where rriiU^^ <^ f^iU^i by canceling the contribution 
of the extra heavy states against that of the extra light ones while keeping the light neutrino 
masses small. This implies a certain level of fine-tuning since extra sterile neutrinos in both 
the heavy and light regimes are necessary and some degree of cancellation between their 
respective contributions is required in order to keep the neutrino masses small. However, in 
such a situation, the contribution of the light extra states could dominate over that of the 
active and induce a rate for the OuPP process larger than the one forecasted from Fig. |2j 
This could thus be a possible solution to an eventual discrepancy between a positive result in 
OuPP decay and a negative result in the searches for neutrino masses in cosmology. Indeed, 
the bounds from cosmology apply to the active SM neutrinos only. 

As an example, we will here consider the Heidelberg-Moscow claim for a positive Oz//3/3 



decay signal 3J]. The accommodation of this signal through only SM neutrinos [see Eq. fl29|) ] 
would require 0.24 eV < ?7i^^ < 0.89 eV at 2a, where the allowed numbers have been ob- 
tained with the ISM results of Sec. Ill Cl following the rather conservative procedure described 
in Ref . 96| . Almost all the error bar comes from the theoretical error of the NME, which 



is much larger than the one associated to the experimental claim. As can be seen in Fig. |21 
the interpretation of this claim as light active SM neutrinos is very disfavoured (see, e.g., 
Ref. 971]) by the constraints from cosmology and neutrino oscillation data. However, this 
signal could be accommodated in a model with heavier neutrinos (which are not bounded 
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by cosmology) mediating the process. Indeed, following Eq. (HT]) . we could reinterpret the 
result as 



0.24 eV < 



heavy 

el 
I 



E ^^^l 



< 0.89 eV. (42) 



This larger contribution to Oz//3/3 decay would not be in conflict with neutrino masses if the 
extra heavy and light neutrino contributions cancel each other in Eq. (HU]) . The level of can- 
cellation required to accommodate the the Heidelberg-Moscow claim with sterile neutrinos 
is actually only at the ~ 50 % level, since the light active neutrinos only fail to explain it by 
about a factor two, given the bounds we have on their masses from cosmology. In Fig. [6] we 
show the degree of cancellation that would be required to accommodate an eventual stronger 



bound on the mass of the lightest neutrino mi from cosmology [35|, |36[ with a discovery of 
0i>l3/3 decay requiring a given m^p when interpreted as the contribution of the SM neutrinos 
alone. The contours for no tension between the mass bound and m^^ are shown together 
with the situations that require cancellations at the level of 50 %, 10 %, 5 % and 1 %. The 
dot represents the present situation of the Heidelberg-Moscow claim that requires a mild 
50 % cancellation to avoid conflict with the cosmological bounds on mi. Notice that, since a 
minimum size of m^^ ~ 10~^ eV is guaranteed for an inverted hierarchy, the required level 
of cancellation never exceeds the 10 % in this situation. On the other hand, if the neutrino 
mass hierarchy is found to be normal, tunings up to ~ 1 % would be necessary to reconcile a 
discovery of m^^ on the same level as the present Heidelberg-Moscow claim with an eventual 
bound of m; < 10^^ eV. 

B. Type-II seesaw models 

In the type-II seesaw, the Standard Model content is expanded by the addition of a scalar 
SU{2) triplet with hypercharge 2 (where the hypercharge is defined such that Q = Y/2+T2): 

I A+/v^ A++ \ 
A= . (43) 

\^ AO -A+/V2J 

The scalar triplet couples to the lepton doublet through the Yukawa terms: 

^ = ifsM - {YA)apL^aiT2ALp + H.C., (44) 

where T2 is the second Pauli matrix. Furthermore, the scalar triplet A has a coupling 
/i to a pair of Higgs fields, so that it gets an induced vev after EW symmetry breaking: 
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FIG. 6: The level of cancellation between the different contributions to the neutrino mass necessary 
to reconcile an eventual discovery of Qu(3(3 decay with size m^^ with a bound on the lightest neutrino 
mass mi from future cosmology or (3 decay experiments. The solid and dashed lines are for normal 
and inverted mass hierarchy respectively. The red, yellow, green, blue and black contours are for 
cancellations to the level of 1 %, 5 %, 10 %, 50 % and no cancellation respectively. The dot 
represents the present Heidelberg-Moscow claim for m^^ and the cosmology bound on mi and 
requires a tuning of ~ 50 % in order to be explained through sterile neutrinos. 



fA = /^w^/2M^, where M^, is the mass of the scalar triplet. The vev of the triplet then 
induces a Majorana neutrino mass 



m^ = 2YaVa = 1a 



Ml 



in Eq. (jH]). This in turn implies that, at low energies, Eq. fl38|) is modified to 

light 



(45) 



(46) 



making m^ the analogous contribution to the one of the heavy neutrinos in the type-I 



seesaw. 



The Oz//3/3 process can be mediated both by the neutrinos and by the scalar triplet (see 
Fig. |3]). However, such charged scalars would have been produced at Tevatron for masses 
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below 100 GeV [98|]. Thus, the contribution of the charged scalar is suppressed with respect 
to the neutrino one. The diagram in the right side of Fig. [3] replaces the neutrino propagator 
by a scalar propagator and thus its amplitude is suppressed by a factor ~ p^ /M\ < 10^^ 
with respect to the SM neutrino contribution. The diagram in the left side of Fig. [3] in 
which one of the W bosons is replaced by the physical charged scalar are also suppressed. 
The scalar is an admixture of the charged components of the Higgs doublet and the scalar 
triplet and its coupling to the quarks is therefore proportional to the quark mass. Thus, the 
amplitude of these contributions turn out to be suppressed by a factor ~ niq/M/^ < 10~^, 
where m^ is the mass of either the up or down quark. Therefore, in this scenario, as in the 
type-I seesaw with all extra states heavy, the light active neutrino contribution dominates 
and the usual description of Oz//3/3 decay in Fig. |5] applies. 

C. Type-III seesaw models 

In the type-III seesaw models the Standard Model is expanded by fermion SU(2) triplets 
with zero hypercharge: 

• (47) 

S- -TP/V2J 

The fermion triplets couple to the SM lepton doublets and the Higgs field through the 
Yukawa terms and have Majorana mass terms of their own: 

=^ = ^SM - ^(Ms).,Tr (S,S^^) - {Y^),J'^%iT2L^ + H.c. . (48) 

The Oz//3/3 decay phenomenology of the type-III seesaw is then completely analogous to 
that of the type-I with the neutral component of the triplet playing the role of the right 
handed neutrino, except that, since the triplet also has charged components, stringent lower 
bounds on its mass exist and in practice only the heavy mass eigenstate regime is available. 
The situation then reduces to the one for the type-II seesaw instead, i.e., the same Oz//3/3 
phenomenology applies, with the replacement 

m^^m^^= ^Y^M-'Y^. (49) 
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D. Mixed seesaw models 

It is interesting to note that the same phenomenology that could stem from a type-I 
seesaw with both heavy and light eigenstates can also arise from a type-11 or III seesaw in 
combination with type-I sterile neutrinos in the light regime. Indeed, adding a type-II or III 
contribution to neutrino masses m^'^ as in Eqs. (145|) and ( 149|) . Eq. ( 132|) would instead read 



771 



A,S 



Ynv/V2 



M.= \ ^ ^ ^ . (50) 

This in turn implies that Eq. ( I36l) is modified to 

light light 

J2 mUl + J2 ^lUli = m^^- (51) 

i I 

Thus, it is possible to have a dominant contribution to Oi//3/3 decay from the extra light 
sterile neutrinos if 7njU^j ^ 7niUl^, while Eq. (ISTl) and the smallness of neutrino masses is 
respected by a cancellation between rrijU'^j and m^'^. The level of the cancellation required 
also corresponds to the one depicted in Fig. O Indeed, the amplitude of Oz//3/3 decay would 
now be 

light light 

Ao^Y. mUlM'"'^^{7ni) + Y, m/f/2,M°^^^(m,) ~ mf,'^M°'^^^(0), (52) 

i I 

and, as an example, the Heidelberg-Moscow claim can be interpreted as 

0.24 eV < \mff\ < 0.89 eV (53) 

in this context. 

V. SUMMARY AND CONCLUSIONS 

We have discussed the general phenomenology of neutrinoless double beta decay (Oz//3/3 
decay) in different types of seesaw models. In particular, we have focused on the contribu- 
tions of the extra degrees of freedom in different mass regimes without assuming preference 
to a particular mass scale. In order to do this, we computed the nuclear matrix element 
(NME) involved in the decay amplitudes as function of the mass of the mediating field, 
detailing all the assumptions performed in each step and estimating the final error due to 
the approximations taken to be at most 30 % for light neutrinos and around 40 % for heavy 
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neutrinos. The results of this computation are pubhcly available at Ref. [48|. In particular, 
the behaviour of the NME is found to be that which can be expected from the propagator 
V(P^ — ?Ti^), where p^ ~ —(100 MeV)^ is the typical momentum transfer between the nucle- 
ons. Thus, the NMEs are essentially constant for m < 100 MeV (light regime) and decrease 
as m~'^ for m > 100 MeV (heavy regime). The transition region around 100 MeV is smooth 
and no significantly new phenomenology takes place at this regime. 

In our discussion we have seen that, for the type-I seesaw, a number of possibilities exist. 
In the case where all the masses of the extra fermion singlets are in the heavy regime, the 
contribution of these states to Oz//3/3 decay is negligible and difficult to constrain directly. 
However, the contribution from the light left handed neutrinos can be rewritten in terms 
of the masses and mixings of the heavy states, which results in quite stringent bounds on 
their mixings. On the contrary, if the extra fermion singlet states are instead in the light 
regime, a GIM-like cancellation with the left handed neutrino contributions occurs. In this 
situation, neutrinos are Majorana particles, but the 0^(3(3 decay rate becomes unobservable, 
being suppressed by at least six orders of magnitude if the mass of the extra states is below 
1 MeV. Thus, considering the extra states alone and neglecting the SM neutrino contribution 
will result in bounds on their mixings which are generally too strong. Such bounds only 
apply when considering extra contributions to neutrino masses beyond the light extra states 
in order to prevent the cancellation from taking place. Dirac neutrinos, where the Majorana 
mass term of the extra fermion states is zero, are a special case of the situation with only 
light extra states. In this scenario, the right and left handed states are exactly degenerate 
and the cancellation is exact. 

If there are extra states in both the light and heavy regimes, then the main contribution 
to the Oz//3/3 transition could come from the light extra states, although some fine-tuning 
is necessary. As such, this could be a way to reconcile a large Oi//3/3 decay rate {e.g., the 
Heidelberg-Moscow claim) with more stringent cosmological bounds with a mild cancellation 
of about 50 %. 

As for the other types of seesaws, current bounds from accelerator experiments place the 
extra degrees of freedom in the heavy regime. This effectively reduces the situation to that 
which appears for the type-I seesaw with only heavy extra states. However, in mixed seesaw 
models, the situation can instead resemble that of the type-I with states in both regimes 
and thus be used to reconcile large 0^(3(3 decay rates with cosmological bounds. 
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In conclusion, the contribution to Oz//3/3 decay from the hght active neutrinos can be 
forecasted by combining present and future neutrino oscillation data on the neutrino mixing 
and mass hierarchy with probes of the absolute neutrino mass scale such as cosmology. 
These predictions can be compared to future Oz//3/3 decay searches so as to gain information 
on the origin and nature of the neutrino masses. In this comparison, we can distinguish the 
following scenarios: 

• The Ouf3/3 process is observed to be in agreement with the forecasted rates. 

This indicates that the light active neutrinos dominate the Oz//3/3 decay rate. Since 
new degrees of freedom are in any event required to give the light neutrinos Majorana 
masses, this implies that there is necessarily new physics above the nuclear scale, so 
that its contribution is suppressed. 

• The OI//3/3 process is observed to be smaller than the forecasted rates. 

This means that there is a partial cancellation between the active and extra neutrino 
contributions. Sterile neutrinos around the nuclear scale are then necessary. A higher 
mass would imply too big a suppression through their NME to show any sizable 
cancellation, while too small masses would make the GIM-like cancellation exact. 

• The Oiyf3(3 process is observed to be larger than the forecasted rates. In this 
situation the light active neutrinos cannot dominate the Oz//3/3 decay rate. Extra sterile 
neutrinos, lighter or around the nuclear scale, could have a significant contribution 
and reconcile the observations. However, the GIM-like cancellation between both 
contributions has to be avoided. This implies either a cancellation between extra 
neutrinos both above and below the nuclear scale (see Sec. IIV A3p or between the 
extra neutrinos and a type-II or III seesaw contribution (see Sec. lIVPp . This is the 
case that would correspond to a confirmation of the Heidelberg-Moscow claim. 

• The 0iy/3(3 process is not observed but was forecasted. While this could imply 
that neutrinos are Dirac and not Majorana particles, it can also be the case that neutri- 
nos are Majorana but extra sterile neutrinos below the nuclear scale are present. Thus, 
the GIM-like cancellation takes place and the Oz//3/3 decay rate becomes unobservable. 

• The Oi//3/3 process is not observed and was not forecasted. This is the most 
pessimistic scenario since it is impossible to draw any conclusion on the nature and 
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origin of neutrino masses. 
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